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O ' Abstract 

^ ■ High-resolution X-ray spectra of SS 433 obtained after a binary egress with the 

Chandra High Energy Transmission Grating Spectrometer (HETGS) were studied. 
Many Doppler-shifted X-ray emission hues from highly ionized elements were de- 



- 

X 



^ i tected. The initial temperature of the jets is estimated to be 20 keV. The lines are 

found to generally be broader than the instrumented resolution. The widths of the Fe 
XXV Ka and Si xill Ka lines correspond to velocity dispersions of 2100 ^340 km s~^ 
and 840 l^gg km s~^ respectively, in terms of Gaussian sigma. Neither the measured 
line widths nor their dependence on the atomic number can be explained by thermal 
broadening alone. Alternative explanations of the observed line widths are discussed, 
including in particular a progressive jet collimation along its axis. 

Key words: Jets — Stars: individual (SS 433) — Stars: binaries: general — 
X-rays: individual (SS 433) — X-rays: spectra 
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1. Introduction 



SS 433 is an enigmatic X-ray binary (orbital period ~13.1 d) with bipolar jets ejecting 
matter at a relativistic velocity of 0.26— times the light speed. The jet axis precesses with 
a period of ~ 162.5 d. Although SS 433 has been studied for more than 20 years since its 
discovery, its fundamental properties, such as the origin of the jet's acceleration and precession, 
and the nature of the compact object, remain unknown. 

The first intensive X-ray studies of SS 433 were performed by the Einstein Satellite. 
Then, SS 433 was revealed to be a variable X-ray source (Grindlay et al. 1984; Band, Grindlay 
1986). Indeed, Watson et al. (1986) and Matsuoka et al. (1986) soon after found a Doppler- 
shifted Fexxv Ka emission line in the X-ray spectra of SS 433 obtained by EXOSAT and the 
Tenma, respectively; they detected one emission line, and Watson et al. (1986) concluded that 
the X-ray emitting part of the receding jet is hidden behind the accretion disk. Thus, the 
length of the X-ray jet was estimated to be shorter than the orbital separation, ~ 10^^ cm. 
From the Ginga data, Kawai et al. (1989) determined the temperature of the jet to be 35 keV 
by modeling the continuum with thermal bremsstrahlung. Brinkmann et al. (1991) applied 
a hydrodynamical numerical model to the Ginga data, and estimated the initial temperature 
of the jet to be as high as ~ 40-70 keV. Since their solution implied a short (~ 10^° cm) jet 
with a high (~ 10^^ cm~'^) density, the receding jet was considered to be easily hidden by the 
accretion disk, to become consistent with the EXOSAT data. 

In 1993, however, a spectrum taken with ASCA, with its high-energy resolution and 
high sensitivity has overthrown the view so far established: the spectrum showed full "pairs" 
of Doppler-shifted emission lines from heavy elements such as Si xiv, S XV, S xvi, Ar xvil, 
Ar XVIII, Ca XIX, Fe xxv, Fe xxvi, and Ni xxvil (Kotani et al. 1994; Kawai 1995). Their 
detection implies that the approaching and receding jets are both visible, so that the X-ray jets 
are longer (~ 10^^ cm) than previously estimated. Kotani et al. (1996) determined the base 
temperature of the jet to be 20 keV from the line flux ratio of Fe xxvi Ka to Fe xxv Ka. 
Further analysis of the spectra revealed that the kinetic luminosity of SS 433 far exceeds the 
Eddington luminosity for any stellar size object (Kotani 1998). Kotani (1998) also suggested 
that the highly ionized iron Ka lines are clearly seen in all of the data sets, while the low-energy 
emission lines from the receding jet are generally weaker than those from the approaching jet, 
and sometimes absent. This can be explained if absorbing matter is present in the line of sight 
to the receding jet. 

Although the ASCA spectra thus revealed many emission lines which no previous mis- 
sion was able to detect, the energy resolution of the ASCA SIS is still insufficient to fully resolve 
the emission line forest; this is particularly the case at energies below ~ 3 keV, where lines 
from different elements overlap and the effect of the absorption is severest. The Chandra High 
Energy Transmission Grating Spectrometer (HETGS) is ideal for an attempt to significantly 
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improve the knowledge obtained with ASCA. An important feature of the HETGS is that spec- 
tra are obtained independently and simultaneously using the High Energy Gratings (HEGs) 
and Medium Energy Gratings (MEGs). The HEGs have a twice-higher spectral resolution and 
a larger effective area above 4 keV than the MEGs, while the MEGs have an extended band- 
pass and a larger effective area for detecting low-energy emission lines. The energy resolution 
(FWHM) of the HEGs is ~ 0.2% at 2.0 keV and ~ 0.5% at 6.0 keV, while that of MEGs 
is ~ 0.14% at 0.8 keV and ~ 0.33% at 2.0 keV. Marshall et al. (2002) reported on the first 
observation of SS 433 using the Chandra/HETGS, made on 1999 September 23. The lines were 
measurably broadened to 1700 km s~^ (FWHM), and the widths did not depend significantly 
on the line-center energy, suggesting that the emission occurs in a freely expanding region of 
constant collimation with an opening angle of 1°.26 ± 0°.06. Furthermore, the excellent angular 
resolution of Chandra has revealed for the first time that the image of SS 433 is extended along 
the east-west direction on a scale of 2"-5", which is comparable to that observed in the radio 
band (Hjellming, Johnston 1981). These impressive results suggest that the Chandra X-ray 
Observatory will give us another surprise regarding this fascinating source. 

2. Observations 

The present observation of SS 433 was performed with the Chandra HETGS on 2001 
May 12 UT10:38-16:33, which was on schedule with our proposal. The total exposure time was 
19.7 ks. The corresponding orbital phase was 0.2832 ± 0.0094 (Gladyshev et al. 1987), when 
the compact object (including the two jets) and the accretion disk emerge out of an eclipse 
by the companion star. The Doppler-shift parameters at this observation were expected to 
be around zto\ ~ +0.075 and z-^a ~ —0.004 for the approaching and receding jets, respectively, 
according to the ephemeris of Margon and Anderson (1989). Usually, the jet directed to the 
east (called "blue jet") is approaching us, and that to the west ("red jet") is receding. However, 
the expected precessional phase on 2001 May 12 was reversed in such a way that the "blue jet" 
is receding and the "red jet" is approaching. This reversal occurs for a minor fraction (~ 50 d) 
of the 162.5 d precession period. 

The observation log is shown in table 1. All the data were acquired in the ACIS-S faint 
mode with CCD chips SO, SI, S2, S3, S4 and S5, and processed in a standard way. The average 
flux of SS 433 in the 1-10 keV band was 8.7 x 10"^^ erg s"^ cm~^, corresponding to a luminosity 
of 2.4 X 10^^ erg s~^, assuming a source distance of 4.85 kpc (Vermeulen et al. 1993). This 
luminosity is typical of SS 433. 
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3. Analysis and Results 



3.1. Data Reduction 

Figure 1 shows the zeroth-order image and the dispersion lines of SS 433 taken with 
the ACIS-CCD. Each of the "arms" of X pattern, formed by diffraction, yields the first-order 
spectrum specified by the grating type (HEG or MEG) and sign of the order (plus or minus). In 
the image, there are no obvious point sources, except for SS 433. On-source event counts have 
been accumulated on the dispersed image over narrow strips along the dispersion lines, of width 
2". 39 on either side. Background counts have been accumulated along similar strips, on both 
sides of each dispersion line, with a distance range of 2". 39-23". 9. The present study utilizes 
the dispersed data of +1 and —1 orders, by co-adding them together because no significant 
discrepancy is seen between them. 

The spectral extraction and data reduction were all performed with the standard pipeline 
for the Chandra HETGS provided by the Chandra X-ray Center (CXO). The following data 
analysis utilizes the Chandra Interactive Analysis of Observations (ClAO) version 2.0 and 
custom routines in XSPEC. The spectra obtained by the HEG/HETGS in the total energy 
band are shown in figure 2. 

When the energy-bin widths are chosen to be comparable to the energy resolution of 
the data, the statistics for each energy bin are so limited that the minimization criterion is 
inappropriate for evaluating the spectral model fits. Instead, the fits are evaluated via minimum 
likelihood method, using the function 

TV 

1=1 

where Xi is the energy at the i-th bin, yi is the observed data at Xj, and y{xi) is the value of the 
fitting function. Minimizing C gives the best-fit model parameters (Cash 1979; XSPEC User's 
Guide for version 11.2.x, p. 151). The C-statistic assumes that the errors on the counts are 
purely Poissonian, and hence it cannot deal with background-subtracted data. Therefore, the 
background is not subtracted in the following analysis. This is justified, because the background 
count rate is only 0.4% of the signal count rate, as is clear from figure 2. 

The jet precession parameters of Margon and Anderson (1989) predict the Fe xxv Ka 
lines from the blue and red jets to appear at energies of 6.236 and 6.718 keV, respectively. 
Indeed, the on-source spectrum exhibits a strong feature over an energy range of 6.3-6.8 keV 
(figure 2). In an expanded spectrum over this energy range, given in figure 3, the emission 
feature is resolved into a series of emission lines. However, none of them corresponds in energy 
to the predicted Fe xxv Ka lines. This suggests that the jet-precession phase has deviated 
from the prediction. Actually, a series of optical spectroscopy of SS 433, performed at Gunma 
Astronomical Observatory (GAO) over a period spanning from 2001 May 12 (simultaneous with 
the present Chandra observation) to the end of 2001 November, consistently imply ~ 20 d phase 
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advance from the prediction (a private communication with K. Kinugasa and H. Kawakita; 
Namiki 2003 in preparation). Below, the jet redshifts are re-evaluated in reference to the 
optical data. 

3.2. Spectral Fitting in the Hard Band 

In figure 3, the HEG/HETGS spectrum in 5.5-8.5 keV energy band around the iron 
lines is fitted with a model consisting of a power-law continuum and narrow Gaussian lines. 
The fitting model is given, as a function of energy E, as 



where c^E) is the photo-electric cross-section at E, is the equivalent hydrogen column 
density, while A and F are the normalization and photon index of the power-law continuum, 
respectively. The fitting model includes Gaussian components corresponding to ten emission 
lines, as specified by the line names in equation (1): four pairs of Doppler-shifted ionized lines 
(Fe XXV Ka, Fe xxvi Ka, Ni xxvil Ka and Fe xxv K/3) from the two jets, and two stationary 
lines (Fe I Ka and Fe I K/3). The rest-frame center energies of the lines were fixed to the 
theoretical values, while the overall Doppler shifts of the blue and red jets were set free; the 
redshift is common among all of the lines belonging to the same jet. The model is subjected to 
photoelectric absorption, with A'^h fixed to the value which is determined in the lower energy 
band (subsection 3.3). 

On the same May 12, the Doppler shifts in the optical band were obtained using Ha 
hues at GAO as -^0.0424 ± 0.0076 and -F0.0204 ± 0.0147 for the blue and red jet, respectively. 
Considering the general agreement between the X-ray and optical Doppler shifts (within < 0.64 
d in the precessional phase; Kotani 1998), the optical redshifts were employed as initial values 
of the X-ray spectral fitting. As a result, the energies of Fe xxv Ka lines were obtained at 
6.410 and 6.524 keV {z = 0.0449 and 0.0267) for the blue and red jets, respectively, as given 
in table 2 and presented in figure 3. These values agree, within the respective errors, with the 
optical measurements at GAO. The observed precessional phase is such that the Fe xxv Ka 
lines from the two jets and the fluorescent iron line from the accretion disk are heavily blended 
together around 6.5 keV. However, thanks to the excellent energy resolution of the HETGS, 
the lines were clearly resolved, except for an overlap between Fe I (6.399 keV) and Fe xxv 
Ka of the blue jet (6.410 keV); these two lines, both known to be generally intense, could not 
be separated even with the energy resolution of the Chandra HETGS. However, this does not 
affect the determination of the jet redshifts, since the other lines contribute. Furthermore, as 
described later, the fluxes of the two overlapping lines can be determined separately with a 
sufficient accuracy by referring to their line widths. 



^-aiE)Nn X X E-^ + Fe I Ka,=o + Fe I K/5,=o 
+ (Fe XXV Ka + Fe xxvi Ka + Ni xxvii Ka + Fe xxv K/?)^ 
-I- (Fe XXV Ka + Fe xxvi Ka + Ni xxvii Ka + Fe xxv K/3) 



^ -^bluG 



(1) 
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The lower panel of figure 3 shows the ratio of the spectrum against the best-fit model 
described by equation (1), with the Gaussians all constrained to be narrow. There, large 
residuals are seen on both sides of each iron emission line. Thus, the lines are inferred to be 
broad, as already reported by Marshall et al. (2002) based on the same instrumentation. As 
shown in figure 4, a much better fit was obtained by allowing the Gaussians for the Doppler- 
shifted lines to have a common finite width; the C-statistics value decreased significantly from 
851.94 (d.o.f. = 595) to 703.57 (d.o.f. = 594), and the residual between the fitting model 
and the data points disappeared. By using Monte Carlo methods to evaluate the absolute fit 
goodness, the narrow Gaussian model was completely rejected, while the broad Gaussian model 
was acceptable at a 2.7 a confidence interval. The Gaussians for the stationary fluorescent lines 
were kept narrow, because they are thought to come from a vicinity of the compact object. 
The width (standard deviation) of the Doppler-shifted lines has been constrained as o" ~ 39.5 
lg;5 eV, or 1850 '^^^ km s~^ in terms of the Doppler velocity dispersion. Although the Fe 
XXV Ka line is intrinsically broad due to the triplet structure (resonance, intercombination and 
forbidden), the Gaussian width decreases insig nificantly to 35.9 tH eV even considering the 
triplet. Here, the intensity ratio of the triplet lines was assumed to be 1.0 : 0.28 : 0.23 for the 
resonance, intercombination and forbidden lines, respectively, calculated by Mewe et al. (1985) 
at a plasma temperature of 20 keV. 

The jet redshifts of SS 433 change continuously with time. Therefore, the line width 
may be caused by drifts in the redshifts during the observation. However, the total data span 
of the present observation, only 21.3 ks, is expected to cause changes in the iron line energies 
by no more than 4.0 eV, which is too small to explain the observed line width. For a further 
confirmation, the spectrum was divided into two parts, the former and latter half of the exposure 
time, and each data were fitted independently. In the case of Fe xxv Ka, the differences of the 
line center energy were obtained as 9.0 ± 29.4 eV and 30.0 ± 41.0 eV for the blue and red jet, 
respectively. As a result, no significant difference was found between the two spectra, and they 
both exhibited broad lines. 

3.3. Spectral Fitting in the Soft Band 

Figure 5 shows the spectra obtained by the MEG and HEG in the soft band (1.5-4.0 
keV). Below, they are fitted jointly. 

First, the data were fitted with a model of a power-law continuum and narrow Gaussian 
lines with photoelectric absorption, just in the same way as was performed in subsection 3.2, 
using the redshifts derived there as initial values. The fitting model is given as 



i-r 



-I- (Si XIII Ka + Si xiv Ka + Si xiv K/3 + S xv Ka + S xvi Ka)^ 
+ (Si XIII Ka + Si xiv Ka + Si xiv K/3 + S xv Ka + S xvi Ka) 



(2) 



6 



where the hehum-hke emission hnes include the resonance and forbidden components, and the 
other annotations of this model are the same as those in subsection 3.2. As a result, a clear 
local minimum in the fit C-statistics has been obtained at zu = 0.0446 and Zrd = 0.0220, both 
close to the input initial values. As shown in figure 5, this fit implies the presence of Si xill 
Ka, Si XIV Ka, S xv Ka, and S xvi Ka lines from the two jet components. 

Although the spectra were thus represented approximately by equation (2), significant 
fit residuals are again observed on both sides of each line. As shown in figure 6, the fit was 
greatly improved by allowing the lines to have a finite common width; the C-statistics value 
has decreased from 1508.82 (d.o.f. = 1236) to 1468.32 (d.o.f. = 1235). The Monte Carlo 
calculation indicates that the former is completely rejected, while the latter is acceptable at 
a 5.0 a confidence interval. The model parameters were obtained as given in table 3. The 
redshift of the blue jet is consistent within the 90% statistical error with that obtained in the 
hard band (table 2). Although the redshift of the red jet is inconsistent with that in the hard 
band, it can be understood by considering the systematic errors (± 0.006 A for HEG; ± 0.011 
A for MEG). The Gaussian width (standard deviation) was obtained to be a = 5.2 t^g eV, or 
870 li5o km s~^ in terms of the Doppler velocity dispersion. 

Even in the broad-line fit in figure 6, there still remain residuals around 1.85 keV, 
approximately exhibiting a double-peaked feature. This feature was successfully modeled by 
two Gaussians with central energies at 1.847 and 1.874 keV and a finite width of cr = 6.8 ^{'^ 
eV. Adding these two lines made the goodness-of-fit of the present model satisfactory at a 1.0 
a confidence interval. However, taking the obtained Doppler shifts into account, the energies 
of these two additional lines do not correspond to the atomic lines of any major element, and 
hence they cannot be represented by the jet components. 

3.4- Spectral Fitting over the Total Band 

By analyzing the HETGS spectra in the Fe Ka and Si Ka bands separately, it has been 
suggested that a common set of 2;bi and Zj-d can reproduce the overall data, but the Fe Ka lines 
are more broadened (in terms of the Doppler velocity dispersion) than the Si Ka lines. In order 
to confirm these two inferences, it is appropriate to fit the two spectral regions simultaneously. 

The data were thus fitted over the total (1.5-9.0 keV) energy band with a model which 
combines those presented in equation (1) and equation (2). The photoelectric absorption, the 
photon index and intensity of the power-law, and the jet redshifts were constrained to be 
common between the two energy bands, while the line widths were allowed to take separate 
values in the hard and soft bands. The unidentified features (at 1.847 and 1.874 keV) were left 
unmodeled in order to reduce the number of free parameters. The results of subsection 3.2 and 
3.3 were used for the initial parameters in the respective energy bands. The joint fit has been 
acceptable (at 5.0 a), and the jet redshifts, line intensities, and the line widths turned out to 
be consistent with those from the individual fits. Although the power-law parameters and the 



7 



absorption column density slightly disagree with those obtained in the limited energy bands, 
this is a relatively minor effect (table 4). 

The most important information obtained from this fit is the relation between the widths 
of low-energy and high-energy lines. As presented in figure 7, the results clearly confirm that 
the velocity dispersion of the high energy lines is significantly larger than that of the low-energy 
lines. 

4. Discussion 

4.I. Jet Temperature 

The intensity ratios among emission lines provide important information on the physical 
conditions of cosmic hot plasmas. In particular, the intensity ratio of the Ka line from Fe xxvi 
to that of Fe xxv provides a clear temperature indicator for very hot plasmas, and hence can 
be used to constrain the jet parameters of SS 433. By comparing the ASCA measurements 
with detailed model calculations, Kotani et al. (1996) estimated the typical initial temperature 
of the jets to be Tq ~ 20 keV. In the present observation, the intensity ratio of the Ka line 
from Fe XXVI to that of Fe XXV was obtained as 0.46 ± 0.17 and 0.42 ± 0.17, for the blue and 
red jet, respectively, which translate to Tq = 19.5lg;g and Tq = 17.2lil'^ keV according to the 
model of Kotani et al. (1996). From these values, the jet activity of SS 433 is inferred to have 
been in a normal state during the observation. 

In contrast, Marshall et al. (2002) derived significantly lower values of Tq, based on the 
Chandra HETGS observation on 1999 September 23 (orbital phase ~ 0.67). In the spectra of 
Marshall et al. (2002), the Fe xxvi lines from the jets were so weak that the intensity ratios of 
the iron lines were 0.30 ± 0.12 and 0.17 ± 0.08 for the blue and red jets, respectively, yielding 
To = 12.9l4;3 To = 8.6 tl'.t keV. Because these temperatures are lower than the usually 
reported values (~ 20 keV), it appears that the high-energy lines, emitted near the base of 
the jet of SS 433, were at that time hidden by something, such as an expanded accretion disk. 
Thus, it is likely that Marshall et al. (2002) observed an unusual low-activity state of SS 433. 
The present paper provides the first report on the high-resolution X-ray spectroscopy of SS 433 
in a relatively normal state, which is close to the past usual states observed with Ginga and 
ASCA. 

4-2. Line Widths 

The Doppler-shifted lines of SS 433, measured with the HETGS, have been successfully 
represented by broad Gaussian models. The obtained line widths are a = 44.5^7^/ eV and a = 
5.0 loi eV, or vpc = 2100 km s"^ and vsi =840 t\lo km s"^ in terms of the Doppler velocity 
dispersion, for the high-energy (mainly iron) and low-energy (mainly silicon) lines, respectively. 
Marshall et al. (2002) also measured finite widths from the Doppler-shifted lines, although the 
widths were not much different from element to element. Figure 8 summarizes the velocity 
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dispersion of various Ka lines, from the present work and Marshall et al. (2002), as a function 
of the atomic number. Thus, the two measurements agree very well on the low-energy lines, 
while not on the high-energy lines. This may be caused by the low activity of SS 433 during the 
observation by Marshall et al. (2002), as suggested by the low values of Tq. If the high-energy 
elements, above calcium, of Marshall et al. (2002) are excluded, figure 8 reveals a clear positive 
dependence of the velocity dispersion on the atomic number. 

In figure 8, the dashed curves represent the velocity dispersion expected when the line 
widths are caused solely by thermal broadening, for several representative ion kinetic tem- 
peratures. Thus, the measured line widths are too large to be explained in terms of thermal 
Doppler effect, for a reasonable range of temperature where the atomic emission lines are sig- 
nificantly emitted. Furthermore, the measured positive correlation of the velocity dispersion 
on the atomic number is opposite to what is predicted by the thermal Doppler effect. For these 
reasons, the thermal broadening is concluded to be inappropriate as an account of the measured 
line widths. 

The Compton scattering is another candidate of generating the line widths. The X-rays 
lines may be narrow when produced in the jets, and then Compton scattered in a surrounding 
medium to get broadened. This may work, because the optical depth of the jets to the electron 
scattering is of the order of unity (Kotani 1998). In this case, a single Compton scattering in 
a medium of electron temperature Tg is expected to shift the X-ray line energy E by AE ~ 
E{4:kTc — E)/mcc'^, where mgC^ is the electron rest-frame energy. Then, the measured values 
of AE ~ 5 eV for silicon {E ~ 1.8 keV) and AE ~ 45 eV for iron {E ~ 6.7 keV) require Te ~ 
0.8 keV and 2.5 keV, respectively. The discrepant values of may be explained by a picture 
that the high-energy lines are produced in a region closer to to the central engine. However, 
such a scattering medium, having the necessary Compton optical depth (~ 1) and the inferred 
low Te, would inevitably produce prominent edges from partially ionized atoms, which are not 
seen in the spectra. Therefore, Compton scattering may not be appropriate, either. 

Another possible origin of the line widths and its dependence on the atomic number is a 
progressive coUimation of the conical opening angle of the jets. Marshall et al. (2002) ascribed 
the line widths to the Doppler broadening that may result from the conical jet outflow, and 
concluded that the emission occurs in a freely expanding region of constant coUimation with an 
opening half-cone angle of O = 1°.23; the opening angle produces a transverse velocity of the 
jet, which causes the obtained widths. Assuming an oppositely directed pair of jets observed 
at an angle a to the line of sight, the Doppler shifts of the blue and red jets are given by 

1 + z = 7{1 ± (3 cos(a)}, 

where P is the velocity of the jet in unit of c, and 7 = (1 — /9^)~^/^. Solving this equation for 7 
as 
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1/2 



Vjet = /3 X C = < 1 



X C, 



where ^av = (-^bi + -2rd)/2, and utilizing the observed redshifts, Wjet = 0.2629 c and a = 92°. 1 are 
obtained. The velocities of the perpendicular component against the direction of jet traveling 
are obtained as v'p^ = sin(a) = 2100 km s~^ and Wgj = vsi sin(a) = 840 km s~^. The velocity 
dispersion can be equated directly with /t'jcttan(6), where / is a form factor depending on 
the emissivity distribution across the jet cross section. Assuming that the density is uniform 
through the cone's cross section and that the component of the velocity which is parallel to the 
jet axis is the same for all fluid elements in the slice, then a simple calculation gives / = 0.74 
(Marshall et al. 2002). Therefore, the present observation yields 



Since the Fe-K lines are thought to be produced close to the central engine while the Si-K lines 
to come from those positions of the jets which are > 10^^ cm away from the center (Kotani 
1998), the present results imply that the jet coUimation is achieved over a distance of > 10^^ 
cm, or > 10^ sec. This interpretation, if correct, provides valuable information as to the jet 
collimation mechanism. 

Since the collimation is thus inferred to occur on a scale much larger than the accretion 
disk, it is more likely to be achieved by magnetic fields, rather than by any funnel-shaped 
structure in the accretion disk. 

4-3. Unidentified Features 

As mentioned in the last paragraph of subsection 3.3, the double-peaked features, which 
can be explained by neither the jet components nor the fluorescent lines from the accretion 
disk, were found around 1.86 keV. This energy corresponds to that of "stationary" Si xill Ka 
line. From a careful look at the spectra in figure 6, similar features are also suggested around 
2.01 keV, 2.46 keV and 2.62 keV, tentatively identified with Si xiv Ka, S XV Ka, and S xvi Ka 
lines. These unidentified features have not been detected for any other observation of SS 433 
using the Chandra HETGS, and their exact nature remains unknown. 

5. Summary 

The analysis of the Chandra HETGS data of SS 433 acquired on the 2001 May 12 have 
yielded the following results: 

1. The observed precessional phase of the jets was shifted by ~ 20 d from the prediction 
based on Margon and Anderson (1989). 
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2. The initial temperature of the blue and red jets are estimated to be 19.5_9;5 and 17.2_g;8 
keV, respectively. 

3. The Doppler-shifted lines were found to be broad, and the corresponding velocity disper- 
sions are 2100 I340 km s~^ and 840ti5o km s~^ for the high-energy and low-energy lines, 
respectively. 

4. While the observed line widths cannot be explained satisfactorily by the thermal Doppler 
effect or the Compton scattering, a successful explanation may be obtained by assuming 
that the line widths originate from the jet opening angle, and that the angle becomes 
narrower as the jets travel away from the central engine. 

M. N. is supported by the Junior Research Associate Program of RIKEN. M. N. also 
would like to thank Dr. H. Matsumoto, Dr. H. L. Marshall, and the members of Center for 
Space Research (MIT) for the study of the Chandra HETGS. We thank the GAO for the good 
optical data and the meaningful discussion and suggestion. 
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Table 1. Observation log. 



Start time (UT) End time (UT) Orbital phase* Exposure time 



2001/05/12 10:38 2001/05/12 16:33 0.2832 ± 0.0094 19.7 ks 



* Calculated from the parameters given by Gladyshev et al. (1987). 



20:00- 




Fig. 1. Zeroth-order image and dispersion lines of SS 433 taken with the ACIS-CCD on 2001 May 12. 
The central point som:ce (SS 433) and 4 lines (grating spectra) are shown in J2000 coordinates. 
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Table 2. Results of the model fit to the HETGS/HEG 
spectrum of SS 433 in the 4.0-9.0 keV band.* 


Parameter 


Blue Red 


Redshift 


0.0460 ± 0.0013 0.0265±0.0014 


Line fluxes^ 
Fe XXV Kq (6.70) 
Fe XXVI Ka (6.97) 
Ni XXVII Ka (7.80) 
Fe XXV K/3 (7.90) 
Fe I Ka (6.40) 
Fe I K/3 (7.06) 


38.5 tit 34.5 tti 

17.6 tl:l 14.3 11? 
9.2 1^;^ < 5.1 

r o +6.2 , 7 +6.1 

6.4 tli 

Q Q +3.3 
■^•O -2.6 


Line width* 


39.5 1^:1 eV 


Power law F 
Power law A 


1.61 ± 0.07 
( 1.63 ± 0.03 ) xlO'^ ph s"^ cm"^ 




1.09 xlO^^ (fixed) 


C-statistic 


703.57 ( d.o.f. = 594 ) 



* All uncertainties refer to statistical 90% confidence limits. 
In the case of redshift, a systematic error of ± 0.0033 must 
be added. 



* The measured line fluxes, in unit of 10^^ ph s^^ cm^^. 
The number in parenthesis shows the rest-frame line energy, 
in unit of keV. 

* Common to all the ion species. 

Column density of photoelectric absorption, determined in 
the soft energy band. 
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Table 3. Results of the model fit to the HETGS / HEG 
+ MEG spectra of SS 433 in the 1.5-4.0 keV band.* 



Parameter 


Blue 


Red 


Redshift 


0.0444 ± 0.0007 0.0220 ± 0.0008 


Line fluxes^ 






Si XIII / Ka 11 (1.85) 


1 +1-2 

I-'-' -1.0 


< 0.9 


Si XIII Ka (1.86) 


n 7 +1.3 
-1.2 


4.2 til 


Si XIV Ka (2.01) 


7 q +1.6 
'••J -1.5 


5.4 til 


Si XIV Kf3 (2.38) 


3.0 tit 


< 1.7 


S XV / Kq II (2.45) 


1-9 tit 


< 1.7 


S XV Ka (2.46) 


o 9 +2.3 
"^•^ -2.0 


9 q +2-2 

-2.0 


S XVI Ka (2.62) 


4.4 til 


2.3 tl^r 


Line width* 




eV 


Power law F 


1.21 ± 


.05 


Power law A 


( 9.41 ± 0.07 ) xlO^ 


^ ph s"^ cm"^ 


Nh 


( 1.09 ± 0.07 ) xlO^^ cm"2 


C-statistic 


1468.32 ( d.o.f. 


= 1235 ) 


* All uncertainties refer to statistical 90% confidence limits. 
In the case of redshift, a systematic error of ± 0.0017 must 



be added. 

'I Forbidden lines. 

The other annotation symbols are the same as in table 2. 
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Table 4. Results of the model fit to the HETGS / HEG 
+ MEG spectra of SS 433 in the 1.5-9.0 keV band.* 



Parameter 


Blue Red 


Redshift 


0.0443 ± 0.0006 0.0227 ± 0.0007 


Line fluxes^ 




Si XIII / Ka 11 (1.85) 


1.0 +H < 0.8 

— i.U 


Si XIII Ka (1.86) 


2.9 tli 4.3 


Si XIV Ka (2.01) 


7.8 t\i 5.8 t\i 


Si XIV Kf3 (2.38) 


3.0 tl.l < 1.5 


S XV / Kq II (2.45) 


2.0 < 1-5 


S XV Kq (2.46) 


3.1 tit 2.8 til 


S XVI Ka (2.62) 


4.5 til 2.4 til 


Line width* 


5.0 eV 

— u.y 


Fe XXV Ka (6.70) 


42.6 +li 30.9 


Fe XXVI Kq (6.97) 


16.9 tli 10.9 til 


Ni XXVII Ka (7.80) 


9.6 til < 5.9 


Fe XXV K/3 (7.90) 


3.4 tli 3.1 tti 


Fe I Ka (6.40) 


6.8 ttl 


Fe I Kf3 (7.06) 


3.4 til 


Line Width" 


44.5 l^^a^ eV 


Power law F 


1.40 ± 0.04 


Power law A 


( 1.18 ± 0.07 ) xlO-2 ph cm-2 


Nu 


( 1.31 ± 0.06 ) xlO^^ cm"2 


C-statistic 


2311.75 ( d.o.f. = 1964 ) 


* All uncertainties refer to statistical 90% confidence limits. 


In the case of redshift, 


a systematic error of ± 0.0017 must 


be added. 




* The width for low-energy lines. 


** The width for high-energy lines. 



The other annotation symbols are the same as table 2. 
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Fig. 2. Chandra HETGS spectra of SS 433 acquired on 2001 May 12, shown for the total (1.0-8.0 keV) 
band. The upper and lower spectra correspond to the source and background, respectively. 
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5.5 6 6.5 7 7.5 8 8.5 



channel energy fkeV) 

Fig. 3. Same on-source spectrum as presented in figure 2, sliown over tlie iron K-line energy band. The 
solid liistograms represent tlie best-fit jet emission model, in which the jet redshifts are left free but the 
lines arc constrained to be narrow. The location of emission lines are labeled with ion species and the jet 
ID; "bl" denotes lines from the "blue jet", and "rd" denotes those from the "red jet". The bottom panel 
shows the intensity ratio against the best-fit model. 




Fig. 4. Same as figure 3 , but the Dopplcr-shiftcd emission lines in the model are allowed to have a 
common width. 
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Fig. 5. Chandra HEG (lower) and MEG (upper) spectra of SS 433 taken on 2001 May 12, shown over the 
sihcon K-hne energy band. The subscripts "r" and "f" are for tlie resonance and forbidden components, 
respectively. The Gaussian lines are constrained to be narrow. The bottom panel shows the ratio against 
the best-fit model. 





\fe I ocity Dispersion of Iron Line (knn/s) 



Fig. 7. Confidence contours for the velocity dispersion of Fe xxv Ka and Si xiii Ka lines. The contours 
correspond to the 99%, 90% and 68% confidence levels from outside to inside. 
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Fig. 8. Atomic number of Ka line vs. its velocity dispersion. The filled and open circles correspond to 
the Chandra HETGS data observed on 1999 September 23 (Marshall et al. 2002) and 2001 May 12 (the 
present work), respectively The dashed lines are velocity dispersions of the Ka lines expected for thermal 
Doppler broadening. The upper one is for an ion kinetic temperature of 70 keV, and lower is for 20 keV. 
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